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Abstract 

A  mathematical  model  for  the  anode  of  a  direct  methanol  fuel  cell  (DMFC)  is  presented.  This  model  considers  the  mass  transport  in  the 
whole  anode  compartment  and  the  proton  exchange  membrane  (PEM),  together  with  the  kinetic  and  ohmic  resistance  effects  through  the 
catalyst  layer.  The  influence  of  key  parameters  on  methanol  crossover  and  anode  performance  is  investigated.  Our  results  indicate  that,  at  low 
current  density  and  high  methanol  concentration,  the  methanol  crossover  poses  a  serious  problem  for  a  DMFC.  The  anodic  overpotential  and 
reaction-rate  distributions  throughout  the  catalyst  layer  are  more  sensitive  to  the  protonic  conductivity  than  to  the  diffusion  coefficient  of 
methanol.  Increasing  the  protonic  conductivity  can  effectively  enhance  the  performance  of  a  DMFC. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  with  proton 
exchange  membrane  (PEM)  as  electrolyte  and  liquid  metha¬ 
nol/water  solution  as  fuel  is  a  promising  power  source  for 
vehicular  applications.  However,  there  are  still  problems 
with  high  anodic  overpotential  and  methanol  crossover  that 
inhibit  its  performance.  Six  electrons  must  be  exchanged  for 
methanol  oxidation,  which  makes  the  oxidation  kinetics 
inherently  slower  [1].  It  has  been  shown  that  Pt-Ru  materials 
possess  substantially  higher  intrinsic  activity  than  platinum 
alone  [2-5],  but  even  with  this  binary  catalyst,  the  anodic 
overpotential  in  a  DMFC  is  still  much  higher  than  that  in  a 
hydrogen  cell.  The  problem  of  methanol  crossover  not  only 
stands  for  the  inefficient  use  of  methanol  but  also  brings 
about  mixed  potential  at  the  cathode  that  adds  further  to  the 
energy  losses. 

Methanol  transport  in  a  DMFC  occurs  by  both  diffusion 
and  convection.  Verbrugge  [6]  developed  a  simple  diffusion 
model  for  methanol  through  a  PEM,  assuming  the  dilute 
solution  theory.  His  model  was  validated  by  experimental 
data,  which  showed  that  the  methanol  diffusion  through  the 
membrane  occurred  nearly  as  readily  as  through  water. 
Cruickshank  and  Scott  [7]  presented  a  simple  model  to 
describe  the  methanol  crossover  in  a  vapor  feed  DMFC 
and  its  effect  on  the  cathodic  overpotential.  The  measured 


*  Corresponding  author.  Tel.:  +886-4-8511217;  fax:  +886-4-8511224. 
E-mail  address:  jeng@mail.dyu.edu.tw  (K.T.  Jeng). 


permeation  rates  of  water  and  methanol  through  a  Nation Ji 
1 17  membrane  under  varied  pressure  differentials  across  the 
PEM  were  used  to  determine  the  essential  parameters  in  the 
model.  This  model  has  also  been  extended  to  include  a  one¬ 
dimensional  model  of  the  potential  distribution  and  concen¬ 
tration  distribution  of  methanol  in  the  anode  catalyst  layer 
for  a  vapor  feed  DMFC  [1,8].  Ren  et  al.  [9]  used  a  DMFC 
to  experimentally  determine  the  electro-osmotic  drag  of 
water  in  the  membrane.  At  higher  current  densities,  a  linear 
relationship  is  found  between  the  cell  current  and  water  flux 
across  the  fully  hydrated  membrane. 

Baxter  et  al.  [10]  presented  a  model  of  the  DMFC  anode, 
which  was  considered  to  be  a  porous  electrode.  The  pores 
were  filled  with  methanol  solution  in  which  all  species  of  the 
reaction  were  allowed  to  transport  in  the  x-direction  through 
the  depth  of  electrode.  Species  movement  in  the  pseudo  y- 
direction  was  taken  into  account  for  water,  methanol,  and 
carbon  dioxide.  This  model  considered  the  catalyst  layer 
only,  however,  a  rather  complex  system  of  equations  was 
necessary  for  this  model.  Scott  et  al.  [11]  described  a  model 
for  the  methanol  transport  processes  that  was  used  to  predict 
the  methanol  concentration  at  the  catalyst  surface  and 
thereby  the  anodic  overpotential.  This  model,  together  with 
an  empirical  model  of  the  open  circuit  voltage  and  the 
cathodic  overpotential  model,  were  used  to  predict  the 
overall  cell  voltage,  current  density  response  of  the  fuel  cell. 

It  is  the  purpose  of  this  study  to  develop  a  mathematical 
framework  to  simulate  the  anode  of  a  DMFC.  This  math¬ 
ematical  model  considers  the  mass  transport  in  the  entire 
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Nomenclature 

Av  specific  area  of  reaction  surface  (cm2/cm3) 

Cch3oh  local  methanol  concentration  (mol/cm3) 

C™H  oh  methanol  concentration  on  anodic  side  of 
PEM  (mol/cm3) 

Qm^oh  methanol  concentration  on  cathodic  side  of 
PEM  (mol/cm3) 

Cph3oh  teecl  methanol  concentration  (mol/cm3) 

Cqj  oh  reference  methanol  concentration  (mol/ cm3) 
Q:l[,0|[  methanol  concentration  at  surface  of  diffusion 
layer  (mol/cm3) 

Coh3oh  methanol  concentration  at  z  =  0  (mol/cm3) 
Ch2o  local  water  concentration  (mol/cm3) 

C™2q  water  concentration  on  anodic  side  of  PEM 
(mol/cm3) 

Q£o  water  concentration  on  cathodic  side  of  PEM 
(mol/cm3) 

CgHjOH  derivative  of  CCh3oh  (mol/cm4) 

Cch^oh  value  of  Cch3oh  evaluated  at  z  =  0  (mol/cm4) 
^ch3oH  effective  diffusion  coefficient  of  methanol  in 
catalyst  layer  (cm2/s) 


Dd 

^ch3oh 

bulk  diffusion  coefficient  of  methanol  in  water 
(cm2/s) 

xyi,eff 

^CH3OH 

effective  diffusion  coefficient  of  methanol  in 
diffusion  layer  (cm2/s) 

^ch3oh 

diffusion  coefficient  of  methanol  in  PEM 
(cm2/s) 

Dm 

^h2o 

diffusion  coefficient  of  water  in  PEM  (cm2/s) 

F 

Faraday’s  constant  (96,500  C/mol) 

i 

protonic  current  density  (A/cm2) 

J*0,ref 

reference  exchange  current  density  (A/cm2) 

i° 

protonic  current  density  at  z  =  0  (A/cm2) 

i 

derivative  of  i  (A/cm3) 

i 

value  of  i  evaluated  at  z  =  0  (A/cm3) 

i 

cell  current  density  (A/cm2) 

kd 

mass  transfer  coefficient  of  methanol  in 
diffusion  layer  (cm/s) 

kf 

mass  transfer  coefficient  of  methanol  from 
feed  stream  to  diffusion  layer  (cm/s) 

km 

mass  transfer  coefficient  of  methanol  in  PEM 
(cm/s) 

Km 

protonic  conductivity  of  ionomer  (S/cm) 

jy-e  ff 

effective  protonic  conductivity  in  catalyst 
layer  (S/cm) 

Ks 

electronic  conductivity  of  solid  phase 
(Pt-Ru/C)  (S/cm) 

Kf 

effective  conductivity  of  solid  phase  in 
catalyst  layer  (S/cm) 

4 

thickness  of  catalyst  layer  (cm) 

4 

thickness  of  diffusion  layer  (cm) 

4n 

thickness  of  PEM  (cm) 

Mh2o 

molecular  weight  of  water  (18  g/mol) 

/Vch3oh 

local  methanol  flux  in  catalyst  layer 
(mol/(cm2  s)) 

^ch3oh 

methanol  flux  through  diffusion  layer 

(mol/(cm2  s)) 

^ch3oh 

methanol  flux  through  PEM  (mol/(cm2  s)) 

/Vdiff 

water  flux  results  from  diffusion  in  PEM 
(mol/(cm2  s)) 

^drag 

water  flux  results  from  electro-osmotic  drag  in 
PEM  (mol/(cm2  s)) 

/Vh2o 

local  water  flux  in  catalyst  layer  (mol/(cm2  s)) 

Nd  ~ 
iVH20 

water  flux  through  diffusion  layer 
(mol/(cm2  s)) 

Nm 

iyn2o 

water  flux  through  PEM  (mol/(cm2  s)) 

Nj 

total  flux  of  water  and  methanol  in  catalyst 
layer  (mol/(cnr  s)) 

R 

universal  gas  constant  (8.3144  J/(mol  K)) 

T 

cell  temperature  (K) 

vd 

superficial  velocity  of  water  in  diffusion  layer 
(cm/s) 

vm 

superficial  velocity  of  water  in  PEM  (cm/s) 

■*CH3OH 

mole  fraction  of  methanol 

Greek  symbols 

aa 

anodic  transfer  coefficient 

y 

order  of  reaction 

£d 

void  fraction  of  diffusion  layer 

ec 

void  fraction  of  catalyst  layer 

4 

volume  fraction  of  solid  phase  (Pt-Ru/C)  in 
catalyst  layer 

4 

volume  fraction  of  ionomer  phase  in  catalyst 
layer 

'/ 

anodic  overpotential  (V) 

electro-osmotic  drag  coefficient  of  water 

Ph2o 

density  of  water  (1  g/cm3) 

<Pm 

potential  of  the  ionomer  phase  (V) 

0s 

potential  of  electronic  conduction  phase  (V) 

anode  compartment  and  the  PEM,  together  with  the  kinetic 
and  ohmic  resistance  effects  through  the  catalyst  layer.  The 
influence  of  the  essential  parameters  on  the  methanol  cross¬ 
over  and  anode  performance  is  investigated. 

2.  Mathematical  model 

Fig.  1  presents  a  schematic  of  the  anode  bonded  to  a  PEM 
in  a  DMFC.  The  system  can  be  divided  into  three  regions:  a 
membrane  region  of  solid  polymer  electrolyte  (PEM),  an 
active  catalyst  region  (catalyst  layer)  that  provides  a  cata¬ 
lytic  site  for  the  oxidation  of  methanol  and  a  diffusion  region 
(diffusion  layer)  composed  of  highly  porous  and  electro¬ 
nically  conductive  material. 

The  diffusion  layer  is  located  adjacent  to  the  flow  channel 
plate  and  its  open  pores  form  transport  avenues  for  the 
distribution  of  reactants  to  the  catalyst  layer  and  collection 
of  carbon  dioxide  into  the  flow  channel.  The  electronically 
conductive  material  conducts  electrons  generated  at  the 
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Fig.  1.  Schematic  diagram  of  the  anode  bonded  to  a  PEM  in  a  DMFC. 


catalyst  layer  into  the  current  collector.  The  catalyst  layer 
is  formed  as  a  thin  film  of  proton  conductive  ionomer 
(e.g.  Nafion  ®)  and  a  carbon-supported  catalyst  (e.g.  Pt-Ru) 
uniformly  dispersed  in  the  ionomer.  The  protons  produced  at 
the  catalyst  layer  are  transported  via  migration  through  the 
PEM.  Moreover,  water  and  methanol  are  also  transported 
through  the  PEM  with  the  protons. 

In  this  study,  the  anode  of  a  single  DMFC  cell  operated 
under  steady-state  isothermal  conditions  is  considered  and  it 
is  assumed  that  the  anode  and  the  cathode  compartments 
have  the  same  pressure.  The  model  equations  are  defined  in 
one  direction  (z-axis)  through  the  anode  and  the  PEM.  The 
origin  is  set  at  the  diffusion  layer/catalyst  layer  interface. 
The  modeling  and  assumptions  for  each  individual  region 
will  be  described  in  the  following  sections. 

2.1.  Methanol  transport  in  the  flow  channel 

The  methanol  concentration  in  the  anode  flow  channel  is 
defined  as  the  feed  concentration.  The  feed  concentration 
decreases  as  the  methanol  solution  flows  down  the  stream. 
As  the  variation  in  the  methanol  concentration  is  small  under 
general  operating  conditions,  one  can  consider  the  feed 
concentration  to  be  constant. 

Carbon  dioxide  is  generated  in  the  anode  electrode  from 
the  electrochemical  reaction.  Gas  bubbles  are  then  formed 
because  of  the  moderate  solubility  of  carbon  dioxide  in  a 
methanol  solution,  and  released  into  the  flow  channel.  The 
released  bubbles  stir  the  flowing  methanol  solution,  which 
complicates  the  methanol  transport  in  the  flow  channel.  The 
methanol  transport  from  the  feed  stream  to  diffusion  layer 
can  be  represented  by 

^CH3OH  =  ^f(^CH3OH  —  ^CH3Oh)  (1) 

where  C^Hj0H  is  the  feed  concentration  of  methanol,  C'(s:H  OH 
the  methanol  concentration  at  the  surface  of  the  diffusion 
layer,  and  k{  is  the  mass  transfer  coefficient.  The  mass 
transfer  at  the  gas-evolving  electrode  was  measured  for 
oxygen  evolution  at  the  vertical  electrode  by  Fouad  and 


Sedahmed  [12].  Scott  et  al.  [11]  correlated  the  oxygen 
evolution  data  and  proposed  that  the  coefficient  of  methanol 
transport  from  the  feed  stream  to  the  diffusion  layer  is 
expressed  using  the  following  empirical  equation 

/  /  \ 032 

kf  =  1.87  x  1CT4  f -  (2) 

\0.003 J 

where  I  is  the  current  density  of  the  fuel  cell.  In  this  study, 
Eq.  (2)  is  adopted  for  the  evaluation  of  methanol  transport 
from  the  flow  channel  to  the  diffusion  layer. 

2.2.  Mass  transport  in  the  diffusion  layer 


The  diffusion  layer  is  composed  of  inactive  carbon  and 
no  electrochemical  reaction  is  expected  within  this  region. 
Therefore,  the  water  and  methanol  flux  remains  constant 
throughout  this  region.  After  passing  through  the  diffusion 
layer,  the  water  is  consumed  partly  by  the  electrochemical 
reaction  in  the  catalyst  layer.  The  remaining  water  then 
migrates  further  through  the  PEM  by  electro-osmosis  and 
diffusion.  The  water  flux  through  diffusion  layer  is  denoted 
as 


Na  =  - 
-,vh2o 


6  F 


■  Nl 


h2o 


(3) 


where  F  is  Faraday’s  constant,  and  N^0  is  the  water  flux 
through  the  PEM,  which  will  be  derived  in  the  PEM 
modeling  in  Section  2.4. 

The  methanol  transport  in  the  diffusion  layer  is  a  combi¬ 
nation  of  diffusion  and  convection;  therefore  the  methanol 
flux  through  the  diffusion  layer  can  be  expressed  as  [13] 


»rd  _  nd.eff  dCCH3OH  .. 

^CH3OH  —  ^CH3OH  +  -XCHsOhNt  (4) 

in  which  Dd’^3fOH  is  the  effective  diffusion  coefficient  of 
methanol  in  the  diffusion  layer,  jich.ou  is  the  local  mole 
fraction  of  methanol,  and  At  is  the  total  water  and  methanol 
flux.  The  effective  diffusion  coefficient  of  methanol  can  be 
derived  from  the  bulk  diffusion  coefficient  DdH3QH  and  the 
void  fraction  of  the  diffusion  layer,  £d,  using  Bruggeman’s 
correction  [14-16]. 


D 


d.eff  _ 
CH3OH  — 


(£d)3/2D 


d 

CH3OH 


(5) 


As  the  concentration  of  methanol  is  low  for  the  dilute 
methanol  solution  used  in  the  DMFC,  the  total  flux  is 
approximately  the  flux  of  water  only,  and  the  local  mole 
fraction  of  methanol  xch3oh  can  be  approximated  by 


C, 


*ch3oh  = 


ch3oh 


c. 


h2o 


CcH: 


OH 


Cch3oh 

Ch2o 


4^h2oCch3oh 

Ph2o 


(6) 


where  MH,o  and  pHl0  are  the  molecular  weight  and  density 
of  water,  respectively.  Eq.  (4)  can  be  rewritten  as 

Ard  _  r,d,eff  dCCH3OH  ,  Mh2oCch3OH  Ard 
^vCH3OH  —  'k'CHsOH  .  I  „  -,vH20  ID 

Ph2o 
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Assuming  ^chjOh  anc*  Ph2 o  t0  constant,  Eq.  (7)  can  be 
solved  within  the  interval  of  <  z  <  0  and  takes  the  form 


/"s  Pi’dAd _ 

„d  LCH3OHe  <-CH3OH,,d  /c,3 

"cH3OH  =  - eVdA«  _  , - ’  (8) 

where  kd  =  OchJoh/ ^d,  is  the  mass  transfer  coefficient  in 
the  diffusion  layer,  Cjl1[  OH  the  methanol  concentration  at  the 
diffusion  layer/catalyst  layer  interface  (i.e.  at  z  =  0),  and 
vd  =  Mh,o^h2o/ Ph,o  is  the  superficial  velocity  of  water  in 
the  diffusion  layer. 

Combining  Eqs.  (1)  and  (8),  one  can  eliminate  Q:Hi0H 
and  express  the  methanol  flux  through  the  diffusion  layer 
using  the  feed  concentration,  C[:Hi0H  as 


n: 


ch3oh 


/"f  ~vd/kd  fO 

uCH3OH  c _ ^CH3OH  d 

evdA’d((vd/A:d)  +  1)  -  1  ' 


(9) 


2.3.  Catalyst  layer 


where  (f>s(z)  is  the  potential  of  the  electronic  conduction 
phase  of  the  catalyst  layer  (i.e.  the  solid  portion  including 
carbon  and  catalyst  particles),  and  cf>m(z)  is  the  potential 
of  the  ionomer  phase.  Both  4>s(z)  and  (f)m  (z)  decrease  in  the 
^-direction,  and  Ohm’s  law  equation  for  each  phase  can  be 
written  as 


d<Ps 

d z 


1 

Kf 


(I~i) 


(15) 


=  _  J_- 
dz  Kf 


(16) 


where  Kf  and  Kf  denote  the  effective  conductivities  of  the 
solid  phase  and  the  ionomer  phase,  respectively.  Therefore, 
the  variation  of  the  overpotential  within  the  catalyst  layer  is 
given  by 


(j_  ,  _LY_  * 

dz  \Kf  +  Kf)1  Kf 


(17) 


The  rate  of  the  electrochemical  reaction  in  the  catalyst 
layer  can  be  described  using  a  Butler- Volmer  rate  expression. 
This  is  then  simplified  to  give  a  Tafel  type  equation  in  terms 
of  the  methanol  concentration  as 


d  i  (  C’ch:!oh\ 

exp 


6aaT;/ 

RT 


(10) 


where  i  is  the  local  protonic  current  density,  A,,  the  specific 
area  of  the  reaction  surface,  «o,ref  the  reference  exchange 
current  density,  y  the  order  of  reaction,  Cch3oh  the  local 
methanol  concentration,  OH  the  reference  methanol 
concentration,  which  is  associated  with  /o,ref,  ota  the  anodic 
transfer  coefficient,  and  t;  is  the  anodic  overpotential. 

In  this  region,  both  water  and  methanol  move  in  the  z- 
direction  due  to  the  diffusion  and  electro-osmotic  drag.  The 
methanol  transport  in  the  catalyst  layer  is  similar  to  that  in 
the  diffusion  layer,  so  the  local  methanol  flux  can  be 
expressed  in  a  similar  fashion  as  Eq.  (7) 


N, 


ch3oh 


„c,eff  dCcH3OH  ,  A7h,oCch3OH 

IS  I 


ch3oh 


di 


P  h2o 


Ni 


h2o 


(11) 


C  eff  r 

where  /JCH  )0h  is  the  effective  diffusion  coefficient  of  metha¬ 
nol  in  the  catalyst  layer,  and  is  the  local  water  flux.  The 
material  balance  for  water  gives 

Nh2o  =  6/7  +  iVH2c  ( 1 2) 


The  flux  of  methanol  in  the  catalyst  layer  decreases  along 
the  z-axis  due  to  its  oxidation  at  the  catalytic  site,  hence  we 
have  the  material  balance  equation  for  methanol  as 
diVcH3OH  _  _  1  di 
dz  6  F  dz 


The  anodic  overpotential  at  any  location  within  the  cat¬ 
alyst  layer  is  defined  as 

n{z)  =  0s(z)  -  (14) 


2.4.  Mass  transport  in  the  PEM 


The  water  transport  through  the  PEM  results  from  both  the 
electro-osmotic  drag  and  diffusion,  i.e. 

N™_0=Ndrdg+Ndiif  (18) 


where  /Vdrag  denotes  the  water  flux  caused  by  the  electro- 
osmotic  drag,  which,  at  a  constant  cell  temperature,  is  propor¬ 
tional  to  the  cell  current  density  /  and  can  be  expressed  as 


A’drag  —  '-H20  .. 


(19) 


where  /.h2o  is  the  electro-osmotic  drag  coefficient  of  the 
water,  which,  for  a  fully  hydrated  PEM,  is  dependent  only  on 
the  temperature.  A^diff  results  from  the  water  concentration 
gradient  in  the  PEM,  and  can  be  expressed  as 


IVdift  =  D 


m 

H20 


/'-'an  /'-'cath 

^h2o  ^h2o 


(20) 


where  £>5t2o  Is  Oic  diffusion  coefficient  of  the  water  in  the 
PEM,  £m  the  thickness  of  the  PEM,  C™20  and  are  the 
water  concentrations  on  the  anodic  and  cathodic  sides, 
respectively.  It  is  appropriate  to  assume  that  the  anodic  side 
of  the  PEM  retains  a  water  saturation  state.  For  a  DMFC 
operated  at  current  densities  higher  than  a  critical  value, 
enhanced  cathodic  water  accumulation  due  to  the  electro- 
osmotic  drag  and  cathodic  reaction  could  raise  the  water 
concentration  on  the  cathodic  side  into  a  fully  hydrated  state 
[9].  For  simplicity,  we  assumed  both  the  anodic  and  cathodic 
sides  of  the  PEM  to  be  fully  hydrated,  eliminating  the  water 
diffusive  transport,  and  A^0  can  be  approximated  by 

=  (2D 


The  methanol  is  transported  through  the  PEM  in  the  same 
way  as  it  is  transported  through  the  diffusion  layer,  so  the 
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methanol  flux  through  the  PEM  (i.e.  methanol  crossover) 
can  be  expressed  in  a  similar  form  as  shown  in  Eq.  (8) 


/'->an  ~vm/km  /'-'cath 

Atm  _  CCH3OH  e _ CCH3OH  m 

iVCH3OH  Qvm/km  _  ] 


(22) 


where  km  =  D™H30H/£m  is  the  mass  transfer  coefficient  of 
methanol  in  the  PEM,  vm  =  Mh2oA^2o/Ph20’  's  t^e  suPer_ 
ficial  velocity  of  the  water  in  the  PEM,  Q'-'h.oh  and  C^OH 
denote  the  methanol  concentration  on  the  anodic  and  catho¬ 
dic  sides  of  the  PEM,  respectively.  Due  to  oxidation  and 
evaporation,  the  methanol  concentration  at  the  cathodic  side 
is  small,  so  we  assumed  that  Cg^0H  is  zero  and  Eq.  (22) 
becomes 


N) 


ch3oh 


/-’an  pvm/km 
uCH3QH  c 

QVm/km  ....  | 


(23) 


3.  Solution  procedure 


By  differentiating  Eq.  (11)  with  respect  to  z  and  combining 
it  with  Eqs.  (12)  and  (13),  we  obtain 
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c.eff  d"CcH3OH  _  Afn2o  (I  —  i 


ch3oh 


d  z2 


N ! 


h2o 


Ph2o  V  6E 

+  |1-*,»Cc»,oh 
Ph2o 


dCcKsOH 

d- 
1  d  i 
6  F  d- 


(24) 


Again,  differentiating  Eq.  (10)  with  respect  to  z  and 
substituting  di]/dz  into  the  so-obtained  equation  with  the 
right-hand  side  of  Eq.  (17)  yields 


d  2i 
dz2 


7  <1Cch3oh  6aaF 
Cch3oh  d-  RT 


1  1  V  1  1 1 di 

K^+Kf)l^Kf\jdz 
(25) 

There  are  two  dependent  variables,  i  and  C’ch3ou-  involved 
in  two  second-order  equations,  Eqs.  (24)  and  (25).  These  two 
equations  can  be  transformed  into  a  fourth-order  system  of 
first-order  equations  as 
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(28) 


(29) 


where  i  and  Cch3oh  stands  for  the  first  derivatives  of  i  and 
Cch3oh,  respectively. 


Four  equations  are  presented  for  the  four  dependent 
variables:  i,  Cch, on.  Cch3oh-  If  their  values  at  z  =  0 

(i.e.  i°,  i  ,  CcHioh’  ancl  CcH,oh)  are  designated,  the  nonlinear 
first-order  initial- value  problem  can  be  solved  to  yield  the  i, 
i,  Cch3oh,  Cch3oh  profiles  using  numerical  methods  (for 
example,  the  Runge-Kutta  method  [17]). 

The  diffusion  layer  is  ionically  insulated,  so  that  the 
protonic  current  density  must  be  zero  at  the  diffusion 
layer/catalyst  layer  interface.  Thus  one  gets 

i°  =  0  (30) 

Employing  the  balance  of  material,  the  methanol  flux  at  the 
diffusion  layer/catalyst  layer  interface  is  equal  to  A^h3oh’ 
and  one  can  write 


Ni 


ch3oh 


=  -D\ 


,eff 


c, 


,0 


ch3oh°ch3oh 


Mh2o  r0  |  N> 

-  Ph  q  Cch3oh^H20  .  6F 


(31) 


Equating  Eqs.  (31)  and  (9),  we  see  that  CCH  OH  is  related 


n  '  "^0 

to  C[' K  ()h|  and  hence  there  are  only  two  initial  values  i  , 
C^h3oh  remaining  to  be  determined. 

At  the  catalyst  layer/PEM  interface,  the  protonic  current 
density  must  be  equal  to  the  cell  current  density  I.  Moreover, 
the  methanol  flux  through  the  diffusion  layer/catalyst  layer 
interface  is  partly  consumed  in  the  catalyst  layer  and  the 
rest  migrates  through  the  PEM  as  methanol  crossover.  The 
chosen  i  ,  C°H;i0H  must  make  the  solution  profiles  satisfy 
the  following  two  conditions 


K  4)  =  I 


(32) 
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2Vch3oh  +  6F  - 
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_iyCH3OHCCH3OH 
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(33) 


in  which  the  methanol  crossover  A^™h3oh 's  evaluated  using 
Eq.  (23).  This  consequently  results  in  a  shooting  problem  in 
terms  of  the  two  initial  values  i  and  Cch3oh>  which  can  be 
solved  using  the  appropriate  shooting  technique.  After 
obtaining  the  solutions  to  the  initial-value  problem,  the 
profiles  of  i  and  Cch3oh  can  be  substituted  into  Eq.  (10) 
to  yield  the  spatial  variation  of  the  overpotential  t]  through 
the  catalyst  layer. 


4.  Results  and  discussion 

Fig.  2  compares  the  calculated  anodic  overpotential  using 
our  model  with  the  experimental  anodic  overpotential  given 
by  Baxter  et  al.  [10].  The  experimental  data  were  obtained 
from  a  methanol  fed,  polymer-electrolyte  fuel  cell  operated 
with  a  hydrogen-producing  cathode  that  was  treated  as  a 
reversible  hydrogen  electrode  (RHE)  by  supplying  it  with 
hydrogen  gas  at  1  atm.  The  methanol  concentration  at 
the  diffusion  layer/catalyst  layer  interface  was  set  to  2  M. 
The  measured  voltages  from  the  methanol/hydrogen  fuel 
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Fig.  2.  Anodic  polarization  of  a  DMFC.  The  experimental  data  correspond 
to  the  methanol  oxidation  data  of  Baxter  et  al.  [10],  Pt-Ru  catalyst 
(90:10  mol%),  2  M  methanol  feed. 

cell  were  corrected  to  reflect  the  anodic  overpotential.  (In 
such  a  fuel  cell,  the  cathodic  overpotential  is  minimized  and 
the  anodic  overpotential  can  be  expressed  as  the  total  over¬ 
potential  subtracted  by  the  ohmic  overpotential  across  the 
membrane.)  Good  agreements  were  obtained  except  at  very 
low  current  densities,  which  is  due  primarily  to  the  Tafel 
kinetic  assumption.  The  parameter  values  used  for  calcula¬ 
tion  are  given  in  Table  1 .  The  effective  values  for  the  catalyst 
layer,  T)C’H  OH,  K™  ,  ,  are  derived  from  their  bulk  values 

using  Bruggeman’s  correction  shown  in  Eq.  (5).  The  rest 
of  them  are  assumption  values  unless  otherwise  noted. 
The  parameter  values  given  in  Table  1  are  regarded  as  the 
base-case  values  in  the  following  simulations. 


0.0  0.2  0.4  0.6  0.8  1.0 
Dimensionless  position 

Fig.  3.  Predicted  variation  of  methanol  concentration  within  the  diffusion 
layer  and  the  catalyst  layer  for  various  current  densities  and  base-case 
conditions. 

Fig.  3  shows  the  methanol  concentration  profiles  in  the 
diffusion  layer  and  the  catalyst  layer  for  various  current 
densities.  As  the  current  density  of  the  cell  increases,  the 
methanol  concentration  decreases  accordingly.  At  lower 
current  densities,  diffusion  dominates  the  mass  transport 
of  methanol  and  only  a  fraction  of  methanol  is  consumed 
in  the  catalyst  layer  by  the  electrochemical  reaction,  therefore 
the  concentration  profiles  are  nearly  linear.  It  should  be  noted 
that,  even  at  high  current  densities,  the  drop  in  the  methanol 
concentration  through  the  catalyst  layer  is  not  as  drastic  as 
the  oxygen  concentration  drop  in  the  cathodic  catalyst  layer 
depicted  in  the  article  by  Bernardi  and  Verbrugge  [18].  This  is 
due  to  the  contribution  of  the  methanol  convection  resulting 


Table  1 

Base-case  parameter  values 


Parameter  Value 


Anodic  transfer  coefficient,  aa  0.14 

Cell  temperature,  T  (K)  353 

Reference  exchange  current  density  times  area,  Avz'o,ref  (A/cm3)  2.5  x  10-4 

Electro-osmotic  drag  coefficient  of  water,  2h2o  3.16  [9] 

Order  of  reaction,  y  0.25 

Catalyst  layer  thickness,  £c  (cm)  0.001 

Diffusion  layer  thickness,  £&  (cm)  0.03 

PEM  thickness,  £m  (cm)  0.015 

Feed  methanol  concentration,  CcH3oh  (mol/cm3)  0.003 

Reference  methanol  concentration,  0H  (mol/cm3)  0.002 

Void  fraction  of  diffusion  layer,  e  0.6 

Diffusion  coefficient  of  methanol  in  water,  D^H30H  (cm2/s)  2.8  x  10-5  [1] 

Diffusion  coefficient  of  methanol  in  PEM,  D™H30H  (cm2/s)  4.9  x  10-6  [1] 

Volume  fraction  of  solid  phase  (Pt-Ru/C)  in  catalyst  layer,  gj;  0.6 

Volume  fraction  of  ionomer  phase  in  catalyst  layer,  0.08 

Void  fraction  of  catalyst  layer,  gc  0.32 

Effective  diffusion  coefficient  of  methanol  in  catalyst  layer,  ^QfeOH  (cm2/s)  5.2  x  10-6 

Protonic  conductivity  of  ionomer,  Km  (S/cm)  0.1416  [1] 

Effective  protonic  conductivity  in  catalyst  layer,  K^f  (S/cm)  3.2  x  10-3 

Electronic  conductivity  of  solid  phase  (Pt-Ru/C),  Ks  (S/cm)  8.13  x  106  [10] 

Effective  conductivity  of  solid  phase  in  catalyst  layer,  Kff  (S/cm)  3.78  x  106 
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Current  density  (A/cm2) 

Fig.  4.  Effect  of  current  density  on  the  methanol  crossover  for  various  feed 
concentrations. 


Fig.  6.  Effect  of  catalyst  layer  thickness  on  the  anodic  overpotential  for 
various  current  densities. 


from  electro-osmosis  with  water,  which  is  absent  in  the 
oxygen  transport  in  the  cathode. 

Increasing  the  current  density  affects  the  methanol  cross¬ 
over  in  two  conflicting  ways.  More  water  migrates  through 
the  PEM,  but  the  methanol  concentration  that  it  takes  along 
decreases,  so  the  methanol  crossover  may  either  increase  or 
decrease  with  the  increasing  current  density.  Fig.  4  illustrates 
the  current  density  effects  on  the  methanol  crossover  for 
various  feed  concentrations.  We  see  that  the  methanol  cross¬ 
over  decreases  with  increasing  current  density  at  feed 
concentrations  lower  than  3  M.  However,  at  4  M  feed  con¬ 
centration,  the  methanol  crossover  increases  slightly  with 
increasing  current  density  then  decrease  after  it  passes 
through  a  maximum. 

If  the  methanol  crossover  is  evaluated  at  a  percentage 
of  the  total  methanol  flux,  we  can  see  the  influence  of  the 
methanol  crossover  on  the  fuel  efficiency  of  a  DMFC.  Fig.  5 
shows  the  methanol  crossover  as  a  percentage  of  the  total 
methanol  flux  transported  from  the  feed  stream  into  the 
electrode  as  a  function  of  the  operating  current  density.  At 


Current  density  (A/cm2) 

Fig.  5.  Methanol  crossover  as  a  percent  of  the  total  methanol  flux  vs.  cell 
current  density  for  various  feed  concentrations. 


lower  current  densities,  a  large  portion  of  the  methanol  fed  into 
the  fuel  cell  is  wasted  because  of  the  methanol  crossover.  For 
feed  concentrations  higher  than  2  M  with  a  current  density 
lower  than  0.1  A/cm2,  the  crossover  is  more  than  60%  of  the 
total  methanol  flux.  The  overall  efficiency  of  a  DMFC  oper¬ 
ated  under  such  a  condition  may  be  no  higher  than  20%. 

Fig.  6  shows  the  effects  of  the  catalyst  layer  thickness  on 
the  anodic  overpotential  for  four  current  densities.  Because 
of  the  increase  in  catalyst  surface  area,  the  overpotential 
decreases  with  increasing  catalyst  layer  thickness  until  a 
certain  thickness,  after  that  it  levels  off  or  even  increases 
(e.g.  0.7  A/cm2)  with  increasing  thickness.  This  is  explained 
by  the  fact  that  a  thicker  catalyst  layer  signifies  a  longer 
distance  through  which  the  current  must  flow  and  the  ohmic 
potential  drop  resulting  in  an  additional  overpotential  across 
the  anode,  especially  when  the  current  density  is  high.  It  is 
seen  that  for  a  DMFC  operated  under  normal  conditions, 
there  may  be  no  advantage  in  having  an  anode  catalyst  layer 
thicker  than  10  pm. 

The  mass  transport  of  the  reactants  and  products  in  the 
catalyst  layer  plays  an  important  role  in  the  electrode 
kinetics  because  it  is  related  to  species  concentration  at 
the  catalytic  site.  The  effective  diffusion  coefficient  of 
methanol  in  the  catalyst  layer,  DC’H  OH,  affects  the  transport 
of  methanol  to  the  catalytic  site  and  the  effective  protonic 
conductivity  in  the  catalyst  layer,  A'JjJ7,  affects  the  proton 
transport  away  from  where  it  is  produced.  Fig.  7  shows  the 
spatial  variations  in  reaction  rate,  di/dz,  throughout  the 
catalyst  layer  at  a  current  density  of  0.3  A/cm2  for  various 
/^ch  Ioh  and  F°r  base-case  value,  the  reaction  rate  is 

higher  in  the  front  portion  (near  z  =  £c)  than  in  the  back 
portion  (near  z  =  0).  The  profile  of  the  reaction  rate  is  only 
slightly  changed  when  Dc|  |  OH  is  reduced  to  one-tenth  of  it  is 
base-case  value,  however,  the  reaction  rate  becomes  fairly 
uniform  throughout  the  catalyst  layer  when  is  increased 
to  10  times  its  base-case  value.  It  is  obvious  that  the  reaction- 
rate  distribution  is  more  sensitive  to  the  protonic  conduc¬ 
tivity  than  to  the  diffusion  coefficient  of  methanol. 
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Dimensionless  position  (z He) 

Fig.  7.  Spatial  variation  of  reaction  rate  throughout  the  catalyst  layer  for 
various  Z)q^0h  and  ■  The  cell  current  density  is  0.3  A/cm2. 

When  the  electrochemical  reaction  occurs  more  uni¬ 
formly  throughout  the  catalyst  layer,  the  catalyst  can  be 
utilized  more  efficiently,  and  this  should  result  in  higher  fuel 
cell  performance.  Figs.  8  and  9  illustrate  the  variation  in  the 
anodic  overpotential  with  the  effective  diffusion  coefficient 
of  methanol  and  protonic  conductivity,  respectively.  From 
the  flat  lines  in  Fig.  8,  we  can  infer  that  the  variation  in  the 
effective  methanol  diffusion  coefficient  in  the  catalyst  layer 
has  almost  no  effect  on  the  anodic  overpotential  of  a  DMFC. 
From  Fig.  9,  however,  we  see  that  the  anodic  overpotential 
can  be  substantially  reduced  if  a  catalyst  layer  with  good 
effective  protonic  conductivity  is  applied.  The  difference  in 
sensitivity  may  result  from  the  fact  that  the  effective  protonic 
conductivity  strongly  governs  the  proton  transport  in  the 
catalyst  layer,  while  methanol  diffusion  is  only  a  portion  in 
methanol  transport.  This  suggests  that,  in  order  to  increase 
the  effective  protonic  conductivity  in  the  catalyst  layer,  the 
void  volume  in  the  catalyst  layer  (space  among  catalyst 
particles)  should  be  filled  with  ionomer  to  a  large  extent. 


Fig.  8.  Variation  of  anodic  overpotential  with  ^cffIoh  f°r  vari°us  current 
densities.  The  circles  indicate  the  base-case  values. 


Effective  protonic  conductivity  KH  (S/cm) 

Fig.  9.  Variation  of  anodic  overpotential  with  AT® f  for  various  current 
densities.  The  circles  indicate  the  base-case  values. 


5.  Conclusions 

We  developed  a  mathematical  model  for  describing  what 
happens  in  the  anode  compartment  associated  with  the 
PEM.  The  formulated  equation  system  was  solved  using 
the  shooting  technique  to  characterize  the  anode  of  a  direct 
methanol  fuel  cell.  The  predicted  polarization  behavior 
compared  well  with  experimental  data  from  previous  work. 
The  results  from  this  study  show  that,  for  a  DMFC  operated 
under  low  current  density  and  high  methanol  feed  concen¬ 
tration  conditions,  a  large  portion  of  the  methanol  fed  into 
the  fuel  cell  is  wasted  due  to  methanol  crossover  making  low 
fuel  efficiency  inevitable.  The  thickness  of  the  catalyst  layer 
affects  the  anodic  overpotential.  Increasing  the  thickness  of 
the  catalyst  layer  may  enhance  the  performance  of  a  DMFC 
to  a  certain  extent,  but  it  is  not  advisable  to  have  a  thickness 
greater  than  10  pm.  The  reaction-rate  distribution  through¬ 
out  catalyst  layer  is  more  sensitive  to  the  effective  protonic 
conductivity  than  to  the  effective  diffusion  coefficient  of  the 
methanol,  hence  increasing  the  effective  protonic  conduc¬ 
tivity  can  effectively  enhance  the  performance  of  a  DMFC. 
This  suggests  that  the  void  volume  in  the  catalyst  layer 
should  be  filled  with  ionomer  to  a  large  extent. 
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